The healing mechanism for excited molecules near metallic surfaces 
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Radiation damage prevents the ability to obtain images from individual molecules. We suggest 
that this problem can be avoided for organic molecules by placing them in close proximity with 
a metallic surface. The molecules will then quickly dissipate any electronic excitation via their 
coupling to the metal surface. They may therefore be observed for a number of elastic scattering 
events that is sufficient to determine their structure. 
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Examining individual molecules at an atomic length 
scale is one of the most challenging modern problems of 
science. An analysis of elastic and inelastic cross sections 
provides a basis for understanding why molecules can 
become fragmented beyond recognition long before the 
time necessary to accumulate sufficient counts for image 
reconstruction. In the case of organic molecules, most 
structural knowledge has been gained using crystalline 
samples, where damaging exposure is distributed over 
many identical molecules, and the image is almost the 
same as that of an undamaged molecule. For molecules 
which cannot be crystallized, there is an urgent need for 
methods that can give information on their structures. 

In single molecule diffraction techniques proposed us- 
ing either femtosecond x-ray pulses [l| or high energy 
electron beams 2|, short wavelengths are used to im- 
age the atomic structure of individual molecules, but 
the radiation doses destroy the molecule after providing 
a diffraction pattern. Imaging with coherent electrons 
from Low Energy Electron Point Sources (LEEPS) is a 
less destructive method to characterize single molecules 
. A holographic interference pattern can be generated 
by LEEPS from which one can reconstruct the object's 
structure. For example, one can reconstruct the structure 
of single polymer strands like DNA ^] . Recently, another 
method using interference patterns with standing x-ray 
waves has been used to map the perylene-tetracarboxylic- 
dianhydride (PTCDA) molecule on a silver surface 0. 
The use of positrons for microscopy is not as well known, 
but it has tremendous potential. Positrons, owing to 
their positive charge, have negative work functions in sev- 
eral materials, and they can be reemitted spontaneously. 
This feature allows the control of slow positron beams 
having energies in the epithermal region. A positron mi- 
croscope for imaging single molecules has been recently 
proposed by Mills and Platzman In order to get 
proper structural information with sufficient resolution. 
Mills and Platzman have estimated that each atom in a 
molecule of 200 atoms should suffer about 50 annihila- 
tions in the course of the positron diffraction. After the 
annihilation, the molecule is in an ionic excited state, 
and under normal circumstances, such a large number of 



ionizations would destroy the molecule. 

Clearly, while exploring the opportunities of imaging 
single molecules with cither x-rays, electrons or positrons, 
it is vital to reduce radiation damage at the level of in- 
dividual molecules. In this letter, we utilize the neutral- 
ization process when a molecule is placed near a metallic 
surface - here called the healing mechanism - and propose 
that it may be applied to heal the radiation damage. The 
main idea is to fill the vacant molecular orbital by an elec- 
tron from the metal in a time short compared to some 
characteristic vibration time of the molecule. Therefore, 
charge neutrality is restored before the onset of the vibra- 
tional modes which would lead to the destruction of the 
molecule. This healing mechanism can be decomposed in 
two main steps. In the first step, when the positive ion is 
sufhciently close to the surface, one conduction electron 
of the metal will tunnel under the influence of the strong 
electric field produced by the ion. In the second step, the 
tunneling electron will fill the hole via an Auger process 
involving an electron-hole pair of the metal. The physics 
of the healing mechanism is closely related to various neu- 
tralization processes [7Lllll [lollTulll.[l^ [l3| and to the 
intermolecular Coulombic decay Il5| . Most of this work 
describes collisions of slow ions with solids. Our study 
involves the neutralization of physisorbed molecules on 
metallic surfaces. Given the wide range of estimates of 
the neutralization time in previous work 8], we must 
perform a conservative estimate to demonstrate the fea- 
sibility of the healing mechanism in our case. 

A schematic of the healing mechanism set-up is shown 
in Fig. n The molecule, in this example PTCDA, is on 
a metallic substrate, which should be chemically inert to 
preserve the integrity of the molecule. Gold, silver or 
other noble metals are possible substrates. In reality, the 
PTCDA geometry might slightly bend near the metallic 
substrate because of the readiness of some PTCDA atoms 
to form weak bonds with the surface |^ . 

When a positive charge is suddenly created on the 
molecule, the resulting electric field attracts the electrons 
of the metal. The probability T that one electron escapes 
from the metal depends on the work function $ and on 
the distance d of the positive charge from the surface. 



FIG. 1: Schematic of a typical healing mechanism set-up: an 
ionized molecule is placed above a metallic substrate and the 
positive charge is at a distance d from the substrate. In this 
example, atoms of a PTCDA fragment are shown. The small 
grey spheres denote H, the large grey spheres denote O and 
the dark large spheres denote C. 



An accurate determination of the tunneling probability 
T would require solving the Schodinger equation, for ex- 
ample using the matching wave function method 
However, the result is mostly dominated by an exponen- 
tial factor given by the WBK approximation ^3 
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dz y/2m{V{z) - EF)/h) , (1) 



where s is the width of the potential barrier and Ep 
is the Fermi energy of the metal. The actual shape of 
the potential barrier V{z) is complicated, but a linear 
approximation of the electrostatic potential at the metal- 
molecule interface yields a simple trian gula r barrier, and 
the penetration probability T becomes 



FIG. 2: Penetration probability T as a function of the dis- 
tance d for a work function "I> = 5 eV. 




FIG. 3: Schematic diagram of the energy levels at the organic 
metal-molecule junction. 



T(d) = exp(~\/2 -$3/2rf2) 
3 



(2) 



where atomic units have been used. Fig. |21 shows that T 
under a barrier characterized by a work function $ = 5 
eV is still about 0.01 for a distance d = 4 A . Actually, 
our major assumption - the triangular barrier - under- 
estimates T since the electric field gets stronger as an 
electron approaches the molecule. Our assumption, as 
we shall see, therefore leads to an overestimate of the 
healing time. 

An electron of the metal is expected to fall into the 
vacant orbital of the molecule. Indeed, a photoemission 
study ^3 of the PTCDA organic molecule on a metal 
surface reveals the energy-level alignment shown in the 



schematic diagram of Fig. where the highest occupied 
molecular orbital (HOMO) level is about 2 eV below Ep. 
A similar diagram of energy levels has been obtained via 
first- principles calculations [20j . 

The energy separation cu (see Fig. ^ between the state 
-04 of the electron spilling out from the metal jl^l and 
the molecular orbital of the hole tpn must be dissipated 
during the recombination. An Auger process involving 
an electron-hole excitation in the metal can mediate the 
capture of the electron in the state ipt into the vacant 
molecular state tjjn- There are several possible electron 
transfer processes between molecules and surfaces jj^ . 
For chemisorbed systems, a hole lifetime is much shorter 
than a typical vibration time |2^ . so that normally the 



hole would be filled before the molecule can fly apart. For 
N2 molecules physisorbed on graphite, experiments sug- 
gest that the neutralization occurs with a characteristic 
time T « lO"-'^'* s ji^. In general, for a molecule far from 
the surface, the corresponding hole may not be filled for 
a while. In this case, the relevant process is the Auger 
neutralization as described by Propst 0. The actual 
calculation of the Auger neutralization rate l/r((i, w) is 
a difficult problem including subtle screening effects at 
the metal- vacuum interface ^8]. Nevertheless, it has been 
suggested by numerical calculations that, as a good ap- 
proximation, the dependence on the distance d and the 
energy uj factorize j2| as 



-(d,a>) = /(rf) — H, 
r tb 



(3) 



where l/rsiio) is obtained using the formula for the 
molecule immersed in the bulk and f{d) is an unknown 
function, usually approximated by an exponential_ 
which decays away from the metal surface. Propst 
however, has shown that the distance of the molecule 
from the surface comes into play mostly through the 
WKB transmission probability. Therefore, we can as- 
sume f{d) « T{d). The matrix element Ma{q) of the 
Auger process in the bulk involves the ground state uq of 
the metal, the excited state Ua of the metal and the den- 
sity operator for excitations with momentum transfer 
q. The actual form of Ma{q) |13 reads 

M„(q) = F{q)Vc{q) < u^\p+\uo >, (4) 

where F(q) is a form factor given by 



^ J '^^ Vr! (r)V't(i")exp(«qr) , 

Vc{q) is the Coulomb interaction Fourier transform 
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and is the volume of the metallic sample. The rate 
of the Auger recombination can be calculated by using 
Fermi's golden rule 



(7) 



where LUao is the energy difference between the metal 
ground-state uq and the excited state Ua ■ Eq. can be 
rewritten in the form 



— = J2S{ci,Lu)V,{qfF{ci)\ 



(8) 



where S'(q, lj) is the dynamical structure factor [2a| de- 
fined by 

5(q,w) = <Uo\pX\uq > \^S{iJ~LJM) ■ (9) 




FIG. 4: Auger rates in the energy range 2 — 10 eV for a metal 
with density parameter rs = 3.1, as in gold. The grey region 
corresponds to ranging from 1 to 00. 



The fluctuation-dissipation theorem relates the dynam- 
ical structure factor S{q, uj) to the dielectric function 
e(q, w) by the formula 
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Substituting Eq.UHlin Eq.|H| we obtain 
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To extract an Auger rate estimate, one can use the Lind- 
hard RPA formula for the dielectric function e{q,Lu). 
Density inhomogeneity effects, due to the lattice and the 
surface, and correlation corrections beyond the RPA are 
not considered in this study. For ipnir), one can use a 
Slater type of orbital (STO) given by 



(12) 



i'nir) = W — exp(-Cr) 

TT 



where 1/C provides the size of the molecular orbital. 
Therefore, assuming that '0t(r) is slowly varying, the 
main form factor contribution to the integral in Ea. llll is 
given by 




5/2 
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(13) 



Fig. El shows typical l/r values in the energy range 2—10 
eV for a gold surface (corresponding to a electron density 
parameter rg — 3.1). Since this rate is proportional to 
the penetration probability T, we show T~^/t in the 
ordinate axis of the plot in Fig 0] The grey region in 
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Fig. 0] corresponds to a STO exponent C ranging from 1 
to oo. For a distance of ~ 4 A, since T is about 0.01 the 
Auger rate becomes of the order of 10^^ 1/s in the lower 
part of the grey region in Fig.^] This value is comparable 
to the highest vibration modes in molecules. All this 
implies that, within distances of less than 4 A, the time 
to fill the hole is estimated to be fast enough to keep the 
molecule from being damaged by molecular vibrations. 
As mentioned previously, the value of the penetration 
probability T given by Eq. Elis a lower bound. This fact 
adds to the robustness of the present model. 

In conclusion, the effects of radiation damage on 
molecules are of concern, because they can be a ma- 
jor obstacle to the emergence of microscopy as a tool to 
determine the structure of delicate biological molecules. 
This is true whether x-rays, electrons or positrons are 
used as the illumination. We have shown that in favor- 
able conditions, ionized molecules on a metal surface can 
quickly dissipate any electronic excitation via their cou- 
pling to the metal. The relaxation is achieved via an 
electron tunneling from the metal to the molecule and 
falling into the vacant orbital. The falling electron gives 
up its energy to an excitation of the metal produced by 
an Auger process. If the molecule is within distances of 
about 4 A from the metal surface, short lifetimes of the 
hole on the molecule (of the order of 10 — 100 fs) prevent 
the destruction of the molecule. Our model may in fact 
overstimate the healing time. This healing mechanism is 
particularly suitable for the positron microscope newly 
proposed by Mills and Platzman 0, but it can be ap- 
plied also to other microscopes imaging single molecules. 
In all these cases, the molecule on the metallic surface 
may be observed without being damaged for a number of 
scattering events sufficient to determine its structure by 
speckle diffraction. These measurements would be ideal 
for studying delicate biological molecules which cannot 
be crystallized. 
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